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ABSTRACT 

We report on the Fermi-LAT observations of the Geminga pulsar, the second 
brightest non-variable GeV source in the 7-ray sky and the first example of a 
radio-quiet 7-ray pulsar. The observations cover one year, from the launch of 
the Fermi satellite through 2009 June 15. A data sample of over 60,000 photons 
enabled us to build a timing solution based solely on 7 rays. Timing analysis 
shows two prominent peaks, separated by A0 = 0.497 ± 0.004 in phase, which 
narrow with increasing energy. Pulsed 7 rays are observed beyond 18 GeV, 
precluding emission below 2.7 stellar radii because of magnetic absorption. The 
phase-averaged spectrum was fitted with a power law with exponential cut-off 
of spectral index T = (1.30 ± 0.01 ± 0.04), cut-off energy Eq = (2.46 ± 0.04 
± 0.17) GeV and an integral photon flux above 0.1 GeV of (4.14 ± 0.02 ± 
0.32) X 10~^ cm~^ s~^. The flrst uncertainties are statistical and the second 
are systematic. The phase-resolved spectroscopy shows a clear evolution of the 
spectral parameters, with the spectral index reaching a minimum value just before 
the leading peak and the cut-off energy having maxima around the peaks. Phase- 
resolved spectroscopy reveals that pulsar emission is present at all rotational 
phases. The spectral shape, broad pulse profile, and maximum photon energy 
favor the outer magnetospheric emission scenarios. 

Subject headings: gamma rays: observations; pulsars: general; pulsars: individual 
(PSR J0633+1746, Geminga) 
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Introduction 
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The Geminga pulsar is the second brightest non-variable GeV 7-ray source in the sky 
and the first representative of a population of radio-quiet 7-ray pulsars . Since its discovery 



as a 7-ray source by SAS-2, more than thirty years ago (JFichtel et al.lll975l : iKniffen et al. 



I975I ). Geminga has been alternatively considered as a unique object or as the prototype of 



a population of hid den dead stars. Fermi has now settled this question with the discovery 
( jAbdo et al.l l2009gl ) of a substantial population of potentially radio-quiet pulsars, of which 
Geminga was indeed the harbinger. 



Geminga was then observed by the COS B 7-ray telescope (JBennett et al.lll977l: 



19811 ). appearing as 2CG 195-1-04 in the seco nd COS B catalog (Swanenburg et al.lll98l[ ) 



and eventually acquiring the name Geminga (JBignami et al.lll983l ). The X -ray source IE 



Masnou et a 



0630 -1-178 detected by the Einstein Observatory in the COS B error box (JBignami et al. 



1983 1) was proposed as a possible cou nterpart, and subseq uently an optical candidate was 
found wi thin the Einstein error box (IBignami et al.lll987l) . which was the bluest object in 



the field (JHalpern fc Tytlerl Il988l : iBignami et al 



1988|) 



The subsequent ROSAT detection of periodic X-rays from this source (JHalpern fc Holt 



9921) prompted a successful search for periodicity in high-energy 7 rays with EGRET 



(IBertsch et al.lll992f ). 
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Geminga has a period of 237 ms and a very stable period derivative of 1.1 x 10^ s s 
that characterize it as a mature pulsar with characteristic age of 3 x 10^ yr and spin-down 
luminosity E = 3.26 x 10^^ erg s~^. 

The determinatio n of the period der i vativ e allowed detectio n of 7-ray pulsation s in the 
previous COS B JBignami fc Caraveol E992[ ) and SAS-2 data jMattox et al.lll992h . Mean- 
while, a high proper motion of 170 mas/yr for the faint my = 25.5 optical counterpart was 
found, confirming the obje ct to be both underluminous and n o more than few hundred pc 
away (IBignami et al.l 119931 ). Using HST, ICaraveo et al.l (Il996[ ) obtained a parallax distance 
for Geminga of 157^34 p c. A comprehensive review of the history of the identification of 
Geminga can be found in IBignami fc Caraveol (119961 ). 



Subsequently, high resolution ast rometry with t 



absolute positioning of Geminga (ICaraveo et al.l 



le Hi pparcos mission allowed for a 40 mas 



19981). Su c h acc urate positioning, together 



with the source proper motion, was used by iMattox et al.l (|l998l ) to improve the quality of 
the timing solution of the pulsar. Recent parallax and proper motion measurements confirm 
the earli er results, yielding a distance of 250lg2° P^ and a proper motion of 178.2 ± 0.4 
mas/yr (JFaherty et al.ll2007l ). 

Analysis of EGRET data showed a double peaked light curve with a peak separation of ~ 0.5 
in phase ( iMayer-Hasselwander et al.lll994J : iFierro et al.lll998l ). The Geminga spectrum mea- 
sured by EGRET was compatible with a power law with a falloff at ~ 2 GeV, but the limited 



-6- 



73 
74 
75 
76 
77 
7B 
79 
80 
81 
82 
83 
84 
85 
86 
87 
88 
89 
90 
91 
92 



EGRET statistics did not allow a nieasurement of the cut-off energy. Deep X-ray observations 
allowed XMM-Newton and Chandra to map the neutron star surface a s it rotates, bringing 



into view different regions contributing differe nt spectral components (ICaraveo et al.l 12004 



De Luca et al.ll2005t iJackson &: Halpern 



trailing the pulsar's motion ( ICaraveo et al. 



20051 



I) as well as an arcmin-sca le bow-shock feature 



20031 : iDe Luca et al.l 120061 ) . A synchrotron ori 



gin of such a non-thermal diffuse X-ray emission trailing the pulsar implies the presence of 
high-energy electrons ( E> 10^^ eV, a value close to the upper energy limit for pulsar wind 
electrons in Geminga) diffusing in a 10 /iG magnetic field. 

Even though Geminga has been one of the most intensively studied isolated neutron stars 
during the last thirty years, it remains of current interest, especially at 7-ray energies where 
its narrow-peaked light curve allows precise timing studies. Thus, it comes as no surprise 
that Ger ainga has been a p rime target for the 7-ray instruments currently in operation: 
AGILE (iTavani et al.l |2009| ) and the Large Area Telescope (LAT) on the Fermi mission 
( JAtwood et al.ll2009l ). Following its launch, the LAT was confirmed to be an excellent instru- 
ment for pulsar studies, observi ng the bright Vela pu lsar ( JAbdo et al.ll2009al ) and discovering 
a variety of new 7- ray pulsars (JAbdo et al.ll2009bl JdJdlJel). including millisecond 7-ray pulsars 
( JAbdo et al.l l2009ir) and a pop ulation of Geminga-like pulsars detected with blind search 
techniques ( lAbdo et al.ll2009g[ l. In this Paper we present the analysis of the Geminga pulsar 
based on the excellent statistics collected during the first year of operations of the Fermi 
mission. 



93 2. 7-ray observations 

94 The Large Area Telescope (LAT) aboard Fermi is an electron-positron pair conversion 

95 telescope sensitive to 7 rays of energies from 20 MeV to > 300 GeV. The LAT is made of 

96 a high-resolution silicon microstrip tracker, a Csl hodoscopic electromagnetic calorimeter 

97 and an Anticoincidence detector for charged particles backgroun d identificatioii . The full 

98 description of the instrument and its performance can be found in lAtwood et al.l (120091 ). 

99 The LAT has a large effective area (peaking at ~8000 cm^ on axis) and thanks to its field of 

100 view (~ 2.4 sr) covers the entire sky every 2 orbits (~ 3 h). The LAT point spread function 

101 (PSF) strongly depends on both the energy and the conversion point in the tracker, but less 

102 on the incidence angle. For 1 GeV normal incidence conversions in the upper section of the 

103 tracker the PSF 68% containment radius is 0.6°. 

104 The data used in this Paper span roughly the first year of operations after the launch of 

105 Fermi on 2008 June 11. The data used for the timing analysis encompass the Launch and 

106 Early Operations (L&EO), covering ~ two months after 2008 June 25, when the LAT was 

107 operated in pointing and scanning mode for check-out and calibration purposes, and extend 
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lOB into the first year of nominal operations up to 2009 June 15. For tlie spectral analysis we 

109 selected only data collected in scanning mode, under nominal configuration, from 2008 Au- 

110 gust 4 to 2009 Jun e 15. We selected ph otons in the 'diffuse' event class (lowest background 

111 contamination, see 



Atwood et al.l l2009l ) and we excluded observations when Geminga was 



112 viewed at zenith angles > 105° where Earth's albedo 7-rays increase the background contam- 

113 ination. We also excluded time intervals when the 15° Region Of Interest (ROI) intersects 

114 the Earth's albedo region. 
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3. Timing Geminga using 7 rays 

Since the end of the EGRET mission, the Geminga timing ephemeris has been main- 
tained using occasional observations with XMM-Newton (Jackson & Halpern 2005; J. Halpern 



priva te communication). While AGILE relied on such X-ray ephemerides ( jPellizzoni et al. 



2OO9I ). LAT densely-sampled, high-precision timing observations yielded an independent tim- 



ing solution. In fact, the LAT timing is derived from a GPS clock on the spacecra ft and times 



of arri val of 7 rays are recorded with an accuracy significantly better than 1 fis f lAbdo et al. 



2009h|). We have constructed a timing solution for Geminga using the Fermi LAT data, ex- 



clusively. For this analysis, we assumed a constant location for the Geminga pulsar calculated 
at the center of the time span of the LAT data set (MJD 54800) using the pos ition reported 



2003). 



bylCaraveo et al.l ( Il998l ) and updated according to the source proper motion ( JFaherty et al. 



We determined an initial, approximate, ephemeris using an epoch-folding search. We then 
measured pulse times-of-arrival (TOAs) by first convertingthe photon event times to a ref- 
erence point at the geocenter using the Fermi sc ience tooljj qtbary, then computing a pulse 
profile using phases generated using TEMP02 (iHobbs et al.l 120061) in its predictive mode. 



Smith et al. 



( 120081 ). This was done with 



The timing accuracy of gtbary was demonstrated in 
~ 22 day segments of data. TOAs were determined from each segment using a Fourier- 
domain cross correlation with a high signal-to- noise template profile. We obtained 16 TOAs 
in this way from 2008 June 25 to 2009 June 15. We fit these TOAs, again using TEMP02, 
to a model with only absolute phase, frequency and frequency first derivative as free param- 
eters. The residuals to the model have an RMS of 251 /is, as shown in Figure H] and the 
model parameters are listed in Table [H The epoch of phase 0.0 given in Table [1] is defined so 
that the phase of the first component of the Fourier transform of the light curve has phase. 
However, in order to assign a smaller phase to the leading peak, we introduced an additional 
phase shift of 0.5 to the timing solution in Table [H Thus, in the light curve shown in Figure 



^http://fermi. gsfc.nasa.gov/ssc/data/analysis/scitools/overview. html 



|2l the epoch of phase 0.0 is the barycentric arrival time MJD(TDB) corresponding to phase 



0.5. 



Table 1: Fermi-L AT Ephemeris for Geminga 



Parameter 



Epoch of position (MJD) 
R.A. (J2000) 
Dec. (J2000) 



Value 



54800 
6:33:54.289 

+17:46:14.38 



Epoch of ephemeris Tq (MJD) 
Range of vahd dates (MJD) 
Frequency f (s^^) 
Freq. derivative / (xlO^^^ s"^) 
Freq. 2nd derivative / (s^'^) 
Epoch of Phase 0.0 (MJD(TDB)) 
Time Units 



54800 

54642 - 54975 

4.21756706493(4) 

-1.95250(9) 



54819.843013078(3) 

TDB 



4. Results 



145 4.1. Light curves 

146 The strong energy dependence of the PSF imposes energy-dependent regions of interest 

147 (ROI) that optimize the signal to noise r atio. Following a p rocedure similar to that used for 
14B the Fermi-LAT pulsar catalogue paper (jAbdo et al.ll2009ll ). to study the pulse profiles we 

149 selected photons within an angle 9 < max[1.6-31ogxo(EGey),l-3] degrees from Geminga. Such 

150 selection provides clean light curves by limiting acceptance of the softer Galactic background. 

151 We used the Fermi tool gtpphase to correct photon a rrival times to the Solar System 

152 barycenter using the JPL DE405 Solar System ephemeris (IStandislj|l998l ) and to assign a 

153 rotational phase to each photon using the timing solution described in Section [31 

154 Figure [2] shows the light curve of Geminga above 0.1 GeV obtained with the energy-dependent 

155 cut. In order to better show the fine structure, we plot the pulse profile using variable- width 

156 phase bins, each one containing 400 events. The photon fiux in each phase interval thus 

157 has a lo" Poisson statistical error of 5%. The dashed line represents the contribution of the 

158 diffuse background, estimated by selecting photons in the phase interval (p = 0.9-1.0 from an 

159 annular region between 2°and 3° from the source rescaled for the solid angle and also taking 

160 into account the energy- dependent selection adopted. The light curve contains 61219 ± 284 
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Fig. 1. — Timing residuals of Geminga derived from the model built using the TOAs of the 
7 rays detected by the LAT (See Section [3] for details). 
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pulsed photons and 9821 ± 99 background photons. 

The pulse profile shows two clear peaks at = 0.141 ± 0.002 (PI) and = 0.638 ± 0.003 
(P2). In order to reveal possible asymmetries in the peaks, we started by fitting the sharp 
peaks with two half-Lorentzian profiles with different widths for the trailing and the leading 
edge. We have chosen this function because it has a simple parameterization and appear 
to fit well the pulse profile of the gamma-ray light curves. We found that Geminga peaks 
show no asymmetries, and PI is broader (FWHM of 0.072 ± 0.002) than P2 (FWHM 0.061 
± 0.001). We also checked if the peaks can be better fitted by a Gaussian profile, finding 
comparable results (PI FWHM of 0.071 ± 0.002) and (PI FWHM of 0.063 ± 0.001), though 
we cannot distinguish between a Lorentz or Gaussian profile. The smallest features in the 
pulse profile appear on a scale of 260 /xs, presumably artifacts of the timing model residuals. 
Figure [2] also contains insets (binned to 0.00125 in phase) centered on the two peaks and on 
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173 the phase interval (j) = 0.9-1.0. This off-peak, or "second interpeak", region contains 789 ± 

174 28 pulsed photons above the estimated background (~ 1.3 x 10^^ of the pulsed flux). This 

175 corresponds to a signal-to-noise ratio of 19cr, indicating that the pulsar emission extends also 

176 in the off-peak, as will be investigated further in Section 4.3. 

Figure [3] shows the pulse profile in five energy ranges (0.1-0.3 GeV, 0.3-1 GeV, 1-3 GeV, 
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Fig. 2. — Geminga light curve above 0.1 GeV using an energy- dependent ROI, shown over 
two pulse periods. The count rate is shown in variable-width bins, each one containing 
400 counts per bin and normalized to 100. Insets show the phase intervals centered on the 
two peaks and on the "second interpeak" region (0 = 0.9-1.0), binned to 0.00125 in phase. 
The dashed line represents the contribution of the diffuse background estimated by selecting 
photons in this "second interpeak" interval in an annulus around the source. 



178 3-10 GeV, > 10 GeV). There is a clear evolution of the light curve shape with energy: 

179 PI becomes weaker with increasing energy, while P2 is still detectable at high energies. 

180 Significant pulsations from P2 are detectable at energies beyond emax ~ 18 GeV, chosen as 
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181 the maximum energy beyond which a x^ periodicity test still attains 6a significance. We 

182 detect 16 photons above 18 GeV, not necessarly coming from the pulsar itself. No particular 

183 features appear at high energies in the bridge region between PI and P2 ("first interpeak"). 

184 Figure S] shows the evolution of the P1/P2 ratio as a function of energy, plotted using 

185 variable-width energy bins. The curve depends very weakly on the bin choice, Figure 4 was 

186 made using 10000 events per bin. A clear decrea sing trend is visib le, as observed in Crab, Vela 



187 and PSR B1951+32 7-r ay pulsars by EGRET ( Thompson! 2004[) and now confirmed for the 



Vela ( JAbdo et al.ll2009al ) and the Crab pulsars (JAbdo et al. 



2010d) by Ferrm LAT. Adopting 



189 the same variable-width energy bins we fit the peaks in each energy range with a Lorentz 

190 function to determine the peak center and width. Figure [5] shows the energy evolution of 

191 the FWHM of PI and P2: both peaks narrow with increasing energy. The decreasing trend 

192 in pulse width of PI and P2 is nearly identical. PI has a FWHM decreasing from 6(f) = 

193 0.098 ± 0.004 to 6(j) = 0.053 ± 0.008, while FWHM of P2 changes from 6^ = 0.092 ± 0.004 

194 to 6(f) = 0.044 ± 0.004 at energies greater than 3 GeV. The decrease in width with energy 

195 does not depend on the shape used to fit the peaks. Figure 8 was made using the Lorenztian 

196 fits, preferred in general because sensitive to asymmetric pulses. While the "first interpeak" 

197 emission is significantly detected up to 10 GeV, emission in the "second interpeak" region 

198 (between 0.9 and 1.0), not detected before, is clearly present at low energies but vanishes 

199 above ~ 2 GeV. 



201 4.2. Energy Spectrum 

202 Spectral analysis was performed using the maximum-likelihood estimator gtlike included 

203 in the standard Fermi Science Tools provided by the FSSC. The fit was performed using a 

204 region of the sky with a radius of 15° around the pulsar position selecting energies between 

205 0.1 and 100 GeV. 

206 We included in the fit a model accounting for the diffuse emission as well as for the nearby 

207 7-ray sources. We modeled the diffuse foreground, including Galactic interstellar emission, 

208 extragalactic 7-ray emission and residual CR background, using the modelcl gll-iem_v02 for 

209 the Galactic part and isotropic_iem_v02 for the isotropic one. 

210 In the fit procedure we fixed the spectral parameters of all the sources between 15° and 20° 

211 from Geminga, and left free the normalization factor of all the sources within 15°. All the 

212 non-pulsar s ources have been m odeled with a power law as reported in the Fermi Bright 

213 Source List (jAbdo et al.ll2009ll ). while all the pulsars have been described by a power law 



^http://fermi. gsfc.nasa.gov/ssc/data/access/lat/BackgroundModels. html 
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214 with exponential c ut-off according to the data reported in the Fermi-LAT pulsar catalogue 



215 flAbdo et al.1 1200911 ). 

216 We integrated the phase- averaged spectrum to obtain the energy ffux. The unbinned gtlike 

217 fit is described by a power law with exponential cut-off in the form: 

^ = iVoE-^ exp (-—) cm-2s-iGeV-i (1) 

dE V ^0/ 

21B where A^'o = (1.189 ± 0.013 ± 0.070) x 10"^ cm^^ g-^ GeV-\ T = (1.30 ± 0.01 ± 0.04) 

219 and Eq = (2.46 ± 0.04 ± 0.17) GeV. The first uncertainties are statistical values for the fit 

220 parameters, while the second ones are systematic uncertainties. Systematics are mainly based 

221 on uncertainties on the LAT effective area derived from the on-orbit estimations, and are of 

222 < 5% near 1 GeV, 10% below 0.1 GeV and 20% above 10 GeV. We therefore propagate these 

223 uncertainties using modified effective areas bracketing the nominal ones (P6_v3_diffuse). 

224 For this fit over the range 0.1 - 100 GeV we obtained an integral photon flux of (4.14 =b 0.02 

225 ± 0.32) X 10"^ cm~^ s~^ and a corresponding energy flux of (4.11 ± 0.02 ± 0.27) x 10~^ 



erg cm ^ s ^. 



We studied alternative spectral shapes beginning with the cut-off function exp[—{E/EQ 



\bl 



The 46 gamma-ray pulsars discussed in lAbdo et al.l (l2010al ) are generally well- described by a 



simple exponential cutoff, 6 = 1, a shape predicted by outer magnetosphere emission models 
(see the Discussion, below). Models where gamma-ray emission occurs closer to the neutron 
star can have sharper "super-exponential" cutoffs, e.g. 6 = 2. Leaving free the exponential 
index b we obtained A'o= (1.59 ± 0.13 ± 0.09) x 10^^ cm-^ g-^ GeV-\ T = (1.18 ± 0.03 
± 0.04), Eo=1.58 ± 0.19 ± 0.1 1) GeV and 6=(0.8 1 ± 0.03 ± 0.06). As previously reported 



for the analysis of Vela pulsar ( lAbdo et al.ll2010bl ). b < 1 can be interpreted by a blend of b 



= 1 spectra with different cutoff energies. Figure [6] shows the results of the phase-averaged 
spectrum in case of b free (dashed line) and b flxed to 1 (solid line). Using the likelihood 
ratio test we found that the hypothesis of b=2 can be excluded since the likelihood of this flt 
being a good representation of the data is much greater than for a power-law flt (logarithm 
of the likelihood ratio being 396) . We have also tried different spectral shapes, like a broken 
power law, but the flt quality does not improve (the logarithm of the likelihood ratio is 212). 



242 4.3. Phase-resolved analysis 

243 We divided the pulse proflle in variable-width phase bins, each one containing 2000 

244 photons according to the energy-dependent cut deflned in Section l4.1[ This choice of bin- 

245 ning provides a reasonable compromise between the number of photons needed to perform 
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246 a spectral fit and the length of the phase intervals, that should be short enought to sample 

247 fine details on the lightcurve, while remaining confortably larger than the rms of the timing 

248 solution (Sec. 3). We have performed a maximum likelihood spectral analysis, similar to 

249 the phase-averaged one, in each phase bin assuming a power law with exponential cut-off 

250 describing the spectral shape. Using the likelihood ratio test we checked that we can reject 

251 the power law at a significance level greater than 5 sigma in each phase interval. Following 

252 the results on phase-averaged analysis of Geminga, we have modeled the spectrum in each 

253 phase interval with a power law with exponential cutoff. Such a model yields a robust fit 

254 with a logarithm of the likelihood ratio greater than 430 in each phase interval. Figure [7] 

255 (below) shows the evolution of the spectral parameters across Geminga's rotational phase. 

256 In particular, the energy cutoff trend provides a good estimate of the high energy emission 

257 variation as a function of the pulsar phase. Table [3] summarizes the results of the spectral 

258 fit in each phase bin. In this case we have fixed all the spectral parameters of all the nearby 

259 7-ray sources and of the two diffuse backgrounds to the values obtained in the phase aver- 

260 aged analysis, rescaled for the phase bin width. 

261 To obtain Fermi-LAT spectral points we divided our sample into logarithmically-spaced 

262 energy bins (4 bins per decade starting from 100 MeV) and then applied the maximum like- 

263 lihood method in each bin. For each energy bin we have used a model with all the nearby 

264 sources as well as Geminga described by power law with fixed spectral index. We have con- 

265 sidered only energy bins in which the source significance was greater than 3a. From the fit 

266 results we then evaluated the integral flux in each energy bin. This method does not take 

267 energy dispersion into account and correlations among the energy bins. To obtain the points 

268 of the Spectral Energy Distributions (SEDs) we multiplied each bin by the mean energy 

269 value of the bin taking into account the spectral function obtained by the overall fit. Figures 

270 [9] to [12] in the Appendix show the SEDs obtained in each phase interval. The fluxes in Y-axis 

271 are not normalized to the phase bin width, whereas in Table [3] of the Appendix the fluxes 

272 are normalized. Figure [7] shows the phase evolution of the spectral index and cut-off energy, 

273 respectively. The spectral index reaches a local minimum around PI (0 ~ 0.14 - 0.15) and, 

274 after a sudden increase, begins to decrease again in the "first interpeak" region, reaching a 

275 minimum of F ~ 1.1 around the leading edge of P2 (0 ~ 0.60 - 0.61). It then starts to rise 

276 again in the phase interval from P2 to the "second interpeak" region (0 = 0.9 - 1.0). 

277 The cut-off energy evolves quite differently as a function of the rotational phase. It closely fol- 



278 lows the pulse profile, thus confirming the observations performed by EGRET (JFierro et al. 



19981 ). which unveiled a correlation between har dness ratio and pulse profile. As shown in 



280 EGRET data and recently confirmed by AGILE (jPellizzoni et al.ll2009l ). the hardest compo- 

281 nent is P2: our phase- resolved scan points to a cut-off around 3 GeV and a spectral index of 

282 ~ 1.0 that become softer through the peak. PI appears to be softer, with a cut-off energy 

283 slightly greater than 2 GeV and a spectral index F ~ 1.2. 
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The phase-resolved spectra show that Geminga's emission in the bridge (or "first int erpeak" ) 
phase in terval ((j) = 0.2 - - 52) is quite differen t from the Crab (jAbdo et al.ll2010d ) or Vela 
pulsars f lFierro et al.l Il998t lAbdo et al.l |2009aJ) . For the Crab pulsar the bridge emission 
shows no evolution and drops to an intensity level comparable to the off pulse emission, 
while for the Vela pulsar it varies substantially but is always seen at high energies. The 
"first interpeak" of Geminga, instead, becomes harder and remains quite strong at high en- 
ergies, as can be also seen in Figure |3l Another difference with respect to the Vela pulsar is 
that G eminga does not ha ve a third peak like the one observed at GeV energies in the Vela 
pulsar jAbdo et al.lboOQal ). 

The analysis of the "second interpeak" region around (f) = 0.9 - 1.0 shows significant emission 
up to ~ 2 GeV (Figure [3]). Moreover the spectrum in this phase interval has been fit with 
a power law with exponential cut-off, obtaining a spectral index F = (1.48 ± 0.17) and Eq 
= (0.87 ± 0.19) GeV, with systematic uncertainties in agreement with those evaluated in 
the phase averaged analysis. A pure power law fit can be rejected with a ~8(T confidence 
level, thus confirming the presence of the cut-off. The presence of the "second interpeak" 
component is also visible in the maps of Figure El where the emission in this phase region is 
not visible at high energies, as expected owing to the spectral cut-off. 
Analyzing the phase evolution of the spectral parameters in Figure [7] it seems that no abrupt 
changes occur in this phase interval and that this emission may be related to the wings of 
the peaks. This fact, together with the newly detected off-peak emission, favors a pulsar 
origin of such "second interpeak" emission, rather than an origin in a surrounding region. 
The detection of off-peak emission, rendered possible by the outstanding Fermi statistics, 
is a novelty of Geminga's high energy behaviour. 



Discussion 



5.1. Light curves and beam geometry 



309 The unprecedented photon statistics collected by Fermi LAT allows for tighter observa- 

310 tional constraints on emission models. The absence of radio emission characterizing Geminga 

311 clearly favors models where the high energy emission occurs in the outer magnetosphere of 

312 the pulsar. 

313 Polar Cap (PC) models, where high energy emission is located near the neutron star surface 

314 ( jDaugherty and Harding II1996I ). are unlikely to explain the Geminga pulsar, since the line 

315 of sight is necessarily close to the magnetic axis for such models where one expects to see 

316 radio emission. 



317 The current evidence against low- altitude emission in 7-ray pulsars (lAbdo et al.ll2009ll ) can 
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31B also be supplemented by constraints on a separate physical origin. In PC models, 7 rays 

319 created near the neutron star surface interact with the high magnetic fields of the pulsar, 

320 producing sharp cut-offs in the few to ~ 10 GeV energy regime. Moreover, the maximum 

321 observed energy of the pulsed photons observed must lie below the 7-B pair production 

322 mechanis m thresh o ld, pr oviding a lower bound to the altitude of the 7-ray emission. Ac- 

323 cording to iBaring] (120041 ). the lower limit for the altitude of the production region r could 

324 be estimated taking advantage of the maximum energy detected for pulsed photons e^ax as 

2 1 

325 r > {emaxBi2/l-76GeV)'' P'tR^^ where P is the spin period, i?* is the stellar radius and B12 

326 is the surface magnetic field in units of lO^^G. For pulsed photons of e^ax ~ 18 GeV, we 

327 obtain rmm > 2.7-R*, a value clearly precluding emission very near the stellar surface, adding 

328 to the advocacy for a slot gap or outer gap ac c eleration locale for the emission in this pulsar. 

329 Outer Gap (OG) models (jCheng et al.lll986l : iRomani Ill996l : IZhang fc Chengll200ll ). where 

330 the high energy emission extends b etween the null charge surface and the light cylinder, 

331 the two-pole caustic (TPC) rn odels ( JDyks and Rudak 1120031 ) associated with slot gap (SG) 

332 ( jMuslimov and Harding 112004 ). where the emission is located along the last open field lines 

333 between the neutron star surface and the light cylinder, or a striped wind model (jPetri 



2OO9I ). where the emission originates outside the light cylinder, could produce the observed 



335 light curve and spectrum. Nevertheless, the observed peak separation of 0.5 is unlikely for a 

336 middle aged pulsar like Geminga in the OG model, if it is true that emission moves to field 

337 lines closer to the magnetic axis as pulsars age. For the OG model this drift leads to < 0.5 
33B peak separations. For TPC models 0.5 peak separation can occur in spite of this shift, that 
339 is, for all ages and spin-down luminosities. 



340 Following the Atlas of 7-ray light curves compiled by IWatters et al.l ( l2009l ). we can use 

341 Geminga's light curve to estimate, for each model, the star's emission parameters, namely 

342 the Earth viewing angle (e with respect to the neutron star spin axis, and the inclination 
angle a between the star's magnetic and rotation axes. Table [2] summarizes the observed 



Model 



a 



Ce 



/i 



n 



TPC 
OG 



30-80,90 
10-25 



90, 55 - 80 

85 



0.7-0.9,0.6-0.8 
0.1-0.15 



Table 2: Earth viewing angles C ki inc lination angles a and beaming factor /q for Geminga, 
as predicted by IWatters et al.l (120091 ) for Outer Gap (OG) and Two Pole Caustics (TPC) 
models. 



344 parameters and gives t he estimated beamin g correction factor /q(«,Cs) 

345 sensitive. It is given by (IWatters et al.ll2009l ): 

/ F^{a; C, <P) sin(C)(iCrf0 



/n(«,Ci 



2jF,{a-XE, 



which is model- 



(2) 
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346 where F^(a;C,0) is the radiated flux as a function of the viewing angle ( and the pulsar 

347 phase (p. In this equation, the numerator is the total emission over the full sky, and the 

348 denominator is the expected phase-averaged flux for the light curve seen from Earth. 

349 The total luminosity radiated by the pulsar is then given by L^ = Att f^FobsD"^ where Fo^s 

350 is the obs erved phase-average d energy flux over 100 MeV and D = 250],]^2o P^ i^ the pulsar 

351 distance (JFaherty et al.l 120071 ) . The estimated averaged luminosity is then L^=3.1x10^''/q 

352 erg s~^, yielding a 7-ray efficiency t]^ = -^ = 0.15fn (rf/lOOpc)^. 

353 Ideally, geometrical values in Table |2] should be compared with independent e stimates, com- 

354 ing e . g. for radio polarization or from the geometry of the pulsar wind nebula (JNg &: Romani 



2004l2008[ ). 



356 Owing to the lack of radio emission, the only geometrical constraints available for Geminga 

357 come from the X-ray observations which have unveiled a faint bow shock structure, due to 

358 the p ulsar motion in the i nterstellar mediura (jCaraveo et al.l 120031 ) and a inner tail struc- 

359 ture (JDe Luca et al.ll2006l : iPavlov et al.l 120061 ). while phase resolved spectros copy yielded a 



360 glimp se of the geometry of the emitting regions as the neutron star rotates (ICaraveo et al. 



2004h . 



362 The shape of the bow shock feature constrains its inclination to be less than 30° with respect 

363 to the plane of the sky. Since such a feature is driven by the neutron star proper motion, 

364 the constraint applies also to the pulsar proper motion vector and th us, presumably, to its 

365 rotation axis, as is the case for the Vela Pulsar (jCaraveo et al.ll200ll ). pointing to an earth 

366 viewing angle ranging from 60 to 90 degrees. 



367 Analysing the pulsar spectral components along its rotational phase, ICaraveo et al.l (120041 ) 

368 concluded that the observed behaviour could be explained in the frame of an almost aligned 

369 rotator seen at high inclination. 

370 However rough, such constraints would definitely favour the OG model pointing to a beam- 

371 ing factor of 0.1-0.15. Such a value turns out to be in agreement also with the heuristic 

372 luminosity law rj ^ ( E'/IO^'^ ) given by lArons I (119961 ) and lWatters et al.l ( 120091 ). that for 



373 the Geminga parameters should yield a value of ~ 17%. For the nominal parallax distance 

374 of 250 pc, a beaming factor of 0.15 would yield a luminosity of L^=4.6xl0'^^ erg s~^. 

375 We note that TPC models, characterized by higher efficiency, would yield higher luminosity 

376 which would account for the entire rotational energy loss for a distance of ~300 pc, well 

377 within the distance uncertainty. On the other hand, a 100% efficiency would translate into a 

378 distance of 730 pc for the OG model, providing a &m limit on the maximum source distance. 
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5.2. Phase resolved spectroscopy 
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The power law with exponential cut-off describes only approximately the phase-averaged 
spectrum of Geminga, since several spectral components contribute at different rotational 
phases. The phase-resolved analysis that we have performed is thus a powerful tool for 
probing the emission of the Geminga pulsar. 

Figure [7] shows a sudden change in the spectral index around each peak maximum. The 
spectrum appears to be very hard in the "first interpeak" region between PI and P2, with 
an index close to F ~ 1.1 and softens quickly after the peak maximum and in the "second 
interpeak" to F ~ 1.5. Caustic models such as OG and TPC predict such behavior as a 
result of the change in emission altitude with energy. Sudden changes in the energy cut-off 
are also predicted, as is also seen for Geminga. Large variations in the spectral index and 
energy cut-off as a f unction of the puls ar phase have already been seen in other p ulsars, such 
as the Crab pulsar jAbdo et al.ll2010ch or PSR J2021-F3651 JAbdo et al.ll2009el ). 
The persistence of an energy cut-off in the " second interpeak" region suggests pulsar emission 
extending over the whole rotation, further supporting the TPC model for Geminga. A similar 
"second interpeak " has been also obse rved by Fermi-LAT in PSR J1836+5925, known as the 
"next Geminga" (JHalpern et al.l 120071 ). Although Geminga is significantly younger, the two 
pulsars share other interesting features, including very similar spectral spectral indexes and 



energy cut-offs in the phase-averaged spectrum, and comparable X-ray spectra ( lAbdo et al. 



2010dl ). 



6. Conclusions 
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In this Paper we presented the analysis of Geminga based on data collected during the 
first year of Fermi operations. The large collecting area of the LAT allows a timing solution 
to be obtained solely from 7-ray data. 

The study of the light curve showed the evolution of the pulse profile with energy, unveiling 
the shrinking of the peaks with increasing energy and providing insights on the highest 
energies with unprecedented detail. Although the phase- averaged spectrum is consistent 
with a power law with exponential cut-off, the phase-resolved analysis showed a much richer 
picture of different spectral components intervening at different rotational phases. The 
phase-resolved analysis has also allowed the detection of the "second interpeak" emission 
indicating a pulsar emission extending over all phases. This feature, never seen before in 
Geminga, was recently also seen by Fermi LAT in PSR J1836+5925 (jAbdo et al.ll2010dl ). 



Our results favor the outer magnetospheric origin for the 7-ray emission. The distance 
uncertainty allows for reasonable values of efficiency for both OG and TPC models, although 
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413 the efficiency for the TPC model becomes too large for distance values just above the nom- 

414 inal one. Future improvements in estimating the distance of Geminga will help to better 

415 strengthen the conclusions and constraining outer magnetospheric models. 

416 The light curve and phase-resolved spectral studies provide a much stronger constraint on 

417 the model geometry. The inclination and viewing angle phase space for peak separation of 
41B 0.5 is very small for the OG, which however provides values compatible with those obtained 

419 from the analysis of Geminga's X-ray behaviour. On the other hand, TPC geometry would 

420 seem more natural for pulsars of Geminga's age that have large gaps. Pulsed emission at all 

421 phases is a common feature of the TPC geometry. It occurs infrequently for OG geometries, 

422 although is present for the large C,e solutions invoked here for Geminga. 

423 The Fermi LAT Collaboration acknowledges generous ongoing support from a number 
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426 Administration and the Department of Energy in the United States, the Commissariat a 
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429 and the Istituto Nazionale di Fisica Nucleare in Italy, the Ministry of Education, Culture, 
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431 (KEK) and Japan Aerospace Exploration Agency (JAXA) in Japan, and the K. A. Wallen- 

432 berg Foundation, the Swedish Research Council and the Swedish National Space Board in 

433 Sweden. Additional support for science analysis during the operations phase is gratefully 

434 acknowledged from the Istituto Nazionale di Astrofisica in Italy and the Centre National 

435 d'Etudes Spatiales in France. 



436 Appendix A: detailed results from phase-resolved spectral analysis 

437 In this Appendix we report all the numerical results and the spectral Energy Distribu- 

438 tions (SEDs) obtained from the phase-resolved spectral analysis of Geminga. Table [3] shows 

439 the spectral parameters obtained from the spectral fit in each phase interval, while Figures 

440 from [9] to [12] show the plots of all the SEDs. 
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Fig. 3. — Geminga light curves in five energy ranges (0.1-0.3 GeV, 0.3-1 GeV, 1-3 GeV, 
3-10 GeV, > 10 GeV). Each light curve is shown over two pulse periods and contains 100 
bins/period. 
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Fig. 4. — Evolution of the ratio P1/P2 with energy, plotted in variable-width energy bins, 
each one containing 10000 events. 
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Fig. 5. — Evolution with energy of the FWHM of PI {bottom) and P2 {top), plotted in 
variable- width energy bins, each one containing 10,000 events. Both peaks narrow at in- 
creasing energies. 
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Fig. 6. — Phase-averaged Spectral Energy Distribution (SED) of the Geminga pulsar. The 
solid line represents the best fit power law with exponential cutoff (i.e b=l), while the dashed 
one represents the best fit power law with exponential cutoff with free exponential index (in 
this case the result is b=0.81). The LAT spectral points (open circles) are obtained using 
the maximum likelihood method described in Section 4.2 
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Fig. 7. — Phase evolution of the spectral index (top) and energy cut-off (bottom) above 0.1 
GeV as the function of the pulse phase, divided in phase bins each containing 2000 photons. 
Vertical bars indicate the combined statistical and systematic uncertainties. For each phase 
interval (defined in Table [3] in the Appendix) a power law with exponential cut-off has been 
assumed. The dashed histogram represents the Fermi-L AT light curve above 0.1 GeV in 
variable-width phase bins of 2000 photons/bin. 
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Fig. 8. — Maps representing the phase interval (0=0.0 - 0.9, top row) compared to the 
second interpeak (0=0.9-1.0, bottom row), in the 2 energy bands 0.1-2 GeV and > 2 GeV. 
Each map represents the photons within 7° from Geminga, binned in pixels of 0.045° {top 
row) and 0.09° {bottom row), smoothed with a gaussian filter with a radius of 2 pixels. In 
the upper left panel we reported the Right Ascension in horizontal axis and the Declination 
in the vertical axis. Bottom row shows that the offpeak point source image is visible at low 
energies but vanishes at E > 2 GeV due to the spectral cut-off. 
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Table 3: Phase interval definitions and corresponding spectral parameters obtained from 
fitting the spectrum with a power law with exponential cut-off. The flux in the third column 
is normalized to the width of the phase bin. The systematic uncertainties are in agreement 
with the ones evaluated for the phase averaged analysis. 
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Fig. 9. — Phase-resolved Spectral Energy Distributions (SEDs) of the Geniinga pulsar in the 
phase range = 0.0 - 0.206. The spectral parameters of each of these spectral distributions 
can be found in Table [31 The fluxes are not normalized to the phase bin width, whereas 
in Table [3] the fluxes are normalized. The curves represent the best fit power law with 
exponential cut-off, while the LAT spectral points (open circles) are obtained using the 
maximum likelihood method described in Section 14.21 
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Fig. 10. — Phase-resolved Spectral Energy Distributions (SEDs) of the Geminga pulsar in 
the phase range (f) = 0.206 - 0.502. 
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Fig. 11. — Phase-resolved Spectral Energy Distributions (SEDs) of the Geminga pulsar in 
the phase range = 0.502 - 0.643. 
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Fig. 12. — Phase-resolved Spectral Energy Distributions (SEDs) of the Geminga pulsar in 
the phase range (p = 0.643 - 1.0. 



